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Abstract 


In this paper, we study the linear complementarity problems on the monotone ex- 
tended second order cones. We demonstrate that the linear complementarity problem 
on the monotone extended second order cone can be converted into a mixed comple- 
mentarity problem on the non-negative orthant. We prove that any point satisfying 
the FB equation is a solution of the converted problem. We also show that the semi- 
smooth Newton method could be used to solve the converted problem, and we also 
provide a numerical example. Finally, we derive the explicit solution of a portfolio 
optimisation problem based on the monotone extended second order cone. 


Keywords: Complementarity problem - Monotone extended second order cone - Portfolio 
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1 Introduction 


The concept of complementarity and complementarity problem, which was firstly intro- 
duced by Karush in [18], is s a cross-cutting area of research and it has a wide range of 
applications in economics, finance and other fields, see [2,3,8,11]. Previous studies show 
that the second order cone programming has played a significant role in complementarity 
problems. The concepts of extended second order cone (ESOC) is introduced by Németh 
and Zhang in [24] and it is a natural extension of the notion of second order cone. Sznajder 
calculated the Lyapunov rank (or bilinearity rank) of ESOC in [27] and proved the irre- 
ducibly of the ESOC. Ferreira and Németh found an efficient numerical method to project 
onto the ESOC [7]. Furthermore, Németh and his collaborators investigated the properties 
of ESOC and used it as a tool for solving various complementarity problems, see [21-25]. 
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They also proposed an application to the optimisation problem of portfolio allocation, 
called the mean-¢? norm (ML2N) model in [28]. The latter paper exhibits advantages 
of the mean-¢? norm (ML2N) model compared to the well-known mean-variance model 
(MV), developed by Markowitz in [19], and the mean-absolute deviation model (MAD), 
introduced in [15]. The application of the ESOC to solving general complementarity prob- 
lems is based on determining its isotone projection sets, concept which is an extension of 
the notion of isotone projection cones (see [14]) and it was introduced in [20]. For the 
importance of the isotone projections in applications see also [13, 26]. 

The importance of the ESOC and ordered vector spaces in investigating and solving 
equilibrium problems important in economics, finance, traffic equilibrium and other fields, 
motivated introducing in in [12] another extension of the second order cone, namely the 
monotone extended second order cone (MESOC). In the latter paper the Laypunov rank 
of MESOC has been determined and it has also been shown that the monotone extended 
second order cone can by used to investigate and solve mixed compelementarity problem. 
Furthermore, Ferreira et. al found a numerical way to project onto MESOC [6] and sug- 
gested applying MESOC to portfolio optimization. In this paper we will show how to solve 
the linear complementarity problem on MESOC and we will give an explicit solution to a 
portfolio optimisation problem on MESOC. 

The structure of the paper is as follows: In Section 2, we introduce the main terminology 
and definitions. In Section 3, we convert the linear complementarity problem on MESOC to 
a mixed complementarity problem on the non-negative orthant. In Sections 4, 5 and 6, we 
will introduce a numerical algorithm which can be used to solve the linear complementarity 
problem on MESOC andin Section 7 we will present a corresponding numerical example. 
Finally, in Section 8 we derive the explicit solution of the considered portfolio optimisation 
problem. 


2 Preliminaries 


Let n > 2 be an integer and R” be the n-dimensional Euclidean space, whose elements are 
identified with column vectors of n components and which is endowed with the classical 
inner product 


(iR xR” > R” 


defined by (x,y) = x'y. Two vectors x, y € R” are called perpendicular if (x,y) = 0, which 
is denoted by x L y. 

If p,q are positive integers such that n = p +q, then for simplicity of notations, we will 
identify the vector space R? x R! with R?+4, by identifying a pair of vectors (x, u) € R? x RY, 
where x € R? and u € R4, with the vector (z',u')' € R?t4. Therefore we will call a pair 
of vector (x, u) shortly vector. Through the above identification the inner product (-,-) in 
R? x IR? becomes 


(x, u), (y, v)) = (x,y) + (u, y), 


for any (x,u), (y, v) € R? x R’. 
In the literature there are various ways of defining cones and various types of cones are 
used. However, in this paper we consider only cones which are closed and convex sets. 


Therefore, for simplicity, we will call a closed set K a cone if and only if ax + By € K, for 
any x,y € K and any a, > 0. A cone K is called proper if it has nonempty interior and 
Kn-K = {0}. 

Let K be a cone. The dual of K is the cone defined by 


K* := {y € R” : (x,y) > 0, Vaz € K}. 
and the complementarity set of K is the set defined by 
CA) := { (x,y): x E€ K, y E K*, x Ly}. 


Definition 1. The monotone extended second order cone (MESOC) is the proper cone 
defined by 


L:= {(x,u) E€ R®x R!’: Tı È T2 È... > Lp = llull} (1) 


Sometimes we will also use the notation £(p, q) to denote that the MESOC is in R? x Rs. 


For the sake of completeness we quote the following four results that will help us proving 
Theorem 5, which are Propositions 3.1, 3.2 in [12] and Propositions 4, 5 in [6]. 


Proposition 1. The dual of the monotone extended second order cone L is the proper cone 
defined by 


j p 
M := {to ER xR’: 5o > 0,Yj € {1,..., p — 1}, X 2 > l} (2) 
i=1 i=1 


From now on, p and q will always denote positive integers, while £ will always denote the 
monotone extended second order cone and M its dual. 


Proposition 2. Let (x,y, u,v) E€ C(L). If u #0,v £0, then 


C(£) = {tema : (x,u) EL, (y,v) EM, 


P Pp 
(x,y) = llull Y to X yi = |/v||,and 3A > 0 such that v = wf 
i=l i=l 


= {tema : (x,u) € L, (y, v) = M, (xi = DD = 0, 
j=l 


p 
Vi =1,...,p— 1, zp = llull, X yi = ||v||,and 3A > 0 such that v = -auh 
i=l 


For any 7 € {1,2...,p}, denote by e’ the vector in R? which has the i-th component one 
and all other components zero and by e the vector in R? with all components one. 


Proposition 3. For arbitrary points (xz,u),(y,v) € R? x R, we have 
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(i) (x,u) € £ if and only if x — |lulle € RS, 
(ii) (y,v) € M if and only if y — |lvlle? € (RELY. 


Proposition 4. Let x,y € R? and u,v € R? \ {0}. Then, we have the following equiva- 
lences: 


(i) (x,0,y,0) E€ C(L) if and only if (x,y) € C(R§4), 
(ii) (w,0,y,v) E€ C(L) if and only if xp = 0, X`} vi = |lvl| and (x,y) € C(RE,), 
0) € C(L) if and only if x; > |lul| for alli, X51 yi = 0 and (x,y) € C(RS,), 


(iv) (t,u,y,v) E€ C(L) af and only if tp = |lull, e) = lull, (ue) = —llullllel], and 


( 
( 
(iii) (x, u,y, 
( 
(x — |lulle, y — llvlle?) € C(RS,). 


Below we list definitions of various types of complementarity problems. 


Definition 2. Let F : R” — R” be an arbitrary mapping and K C R” an arbitrary cone. 
The complementarity problem defined by K and F is 


find an x E K, such that 


CUE ae € C(K) 


If T € R"*” is a constant matriz, r € R” is a constant vector and F(x) = Tx +r, then 
the problem CP(F,K) is called the linear complementarity problem defined by T, r, and K 
and it is denoted by LCP(T,r, K). 


Definition 3. Let G: R? x R! > R’, H : R? x R45 RY, F : R? x R! — R? be arbitrary 
mappings and K C R? an arbitrary cone. The mized implicit complementarity problem 
defined by K, G, H and F is 


find an (x,u) € R? x R! such that 
H(x,u) = 0, (F(x,u), G(x, u)) € C(K) 


Definition 4. Let G: R? x R! > R’, H : R? x R! > R! be arbitrary mappings and K an 
arbitrary cone. Then, the mixed complementarity problem defined by G, H and K is 


MiICP(G, H, F, K) := i 


find an (x,u) E R? x R1, such that 
H(x,u) = 0 and (x,G(z,u)) € C(K) 


MiCP(G, H,K) := i 


3 The Linear Complementarity Problem on the MESOC 


Theorem 5. Let (x,u), (y, v) be arbitrary vectors with x,y € R? and u,v € R1. Consider 


the nonsingular block matrix 
T= A B 
=- (C D)’ 


where A € R?*?, B € R°*4, C e RI*? and D € RI% are constant matrices. Then, for 
arbitrary vectors z* and r, such that z = (x,u) and r = (y,v), the following statements 
hold: 


(i) 


ia 
2. 
Ss 

© 


(v) 


Let u = 0. Then, z is a solution of LCP(T,r,£L) if and only if x is a solution of 
LCP(A,y, R£) , tp = 0 and X$ (Az: + yi) = ||Cx + vll. 


Let Cx + Du +v = 0. Then, z is a solution of LCP(T,r, £L) if and only if x is a 
solution of MiCP(G,H,R&,) , ti > |lull, and X4- (4r + Bu + v); = 0, where G 
and H are defined by the formulas G(x', u’) = Ax' + Bu' +y and H(x',u') = 0. 


Letu #0 # Cx + Du+v. Then, z is a solution of LCP(T,r, L) is equivalent to z is 
a solution of MiỌICP(G, H, F, RE), where F, G and H are defined by the formulas 


F(x u) =a! — |lu'|le, G(x, u) = Av’ + Bu’ +y—||Ca' + Du’ + vlje, 


and 
H(x' u) = u'e! (Ax' + Bul +y) + |lu'||(C2' + Du’ +). 


Denote z = (Z,u) = (x — ||u|je,u) and let u #0 A Cr+ Du+v. Then, z is a solution 
of LCP(T,r,£) is equivalent to Z is a solution of MiCP(G,H,R&_,), where G and 
H are defined by the formulas 


G(x, u) = A(x + |lu'le) + Bul + y — ||Ca’ + |lu'|le) + Du’ + vlje, 
and 
H(x' u") = u'e! (A(x + |lu’|le) + Bul +y) + |Ju’||(C(a’ + |Ju'lle) + Du’ + v). 
When u #404 Cx + Du +v, the problem of finding a solution z = (x,u) of the linear 


complementarity problem LCP(T,r,£) is converted to a problem of finding a vector 
z = (x,u) such that (a, 8) € C(R’.), where 


T1 — T2 (Ax + Bu +y) 
rg — T3 52 (Ar + But y)i 
a= : and B = : 
Tp- — Tp P-i(Ax+ Bu +y); 
Lp — llull im(4r + Bu +y); 


Moreover, denote 
p—1 p—1 
ziw) =X w +e, = X w + llul, 
j=i j=l 


for any i = 1,2,...,p— 1 and any z',w' € RP, w € R°. Let x,(w') = |lu'||. Then, 
the problem of finding a vector z = (x,u) such that (a, 6) € C(R") is equivalent to 
the problem of finding a solution of MiCP(Ĝ, H, RY 1), where 


(Az'(w’) + Bu’ +y) 
yn | Da (Ae! (w) + Bul +y): 


Di (Aa'(w') + Bul +y): 
and i 
H (w, u) = u'e! (Ax' (w) + Bul + y) + llul (Cx' (w) + Du’ +v) 
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(vi) Let t = |lul|. Then, z is a solution of LCP(T,r,£) if and only if x is a solution of 
MiCP(G, H,R*’), where G and H are defined by the formulas 


(Ax'(w, t) + Bu +y) 


2 
pad s Ag’ w’, A T Bu’ + i 
G(w',u’, t’) = Vigil ( l ) y) é Re L 


PTT (Aa! w, t) + Bu +y); 


7 al Ifat H zl 1 Vlt F! Du! 
uad- u'e (AEB TY E a U u +v) 
t? — |[u'| 
and 
wi twt... +w Ht 
Wwy Foa HU FE 
x'(w, t) = : 
wp HË 
t 
Proof. 


(i) By the definition of the linear complementarity problem, z = (x,0) is a solution 
of LCP(T,r,£) if and only if (7,0, Ar + y,Cx+v) € C(L), which, by using item (ii) in 
Proposition 4, is equivalent to x, = 0, \{}_\(Azit+y) = ||Cr+vl| and (x, Arty) € C(R§,). 
Finally, that is further equivalent to x being a solution of LCP(A,y, R&,). 

(ii) Let Cr + Du +v = 0. By the definition of the linear complementarity problem, 
z = (x,u) is a solution of LCP(T,r, £) if and only if (x, u, Az + But y,0) € C(L), which, 
by using item (iii) of Proposition 4, is equivalent to x; > |lul|, e'(Av + Bu + y) = 0 and 
(x, Ax + Bu +y) € C(R&,). We conclude that z = (x,u) is a solution of LCP(T,r, £) if 
and only if z = (x,u) is a solution of MiCP(G, H,R&,). 

(iii) By using the definition of linear complementarity problem, if z = (x, u) is a solution 
of LCP(T,r, £), then we have (x,u, Ax + Bu + y,Cx + Du +v) € C(L). Then, from 
item (iv) of Proposition 4 and the equality case of the Cauchy inequality, we have that 
(x,u, Ax + Bu + y, Cx + Du +v) € C(L) is equivalent to the existence of a A > 0 such 
that the following equations hold: 


£p = lull, 
Cr + Du +v = —\u, (3) 
e' (Az + Bu +y) = ||Cz + Du + v|| = Alļjul| (4) 
and 
(x — |lulle, Ax + But y— ||Cx + Du + ve?) € CRZ4). (5) 


By using (5), we conclude that 
(F(x, u), G(x,u)) € C(R§,). 


By using equation (3) and (4), we have 


H(az,u) = ue! (Az + Bu +y) + |lul|(Cz+ Du +v) =0. 
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Thus, z being a solution of LCP(T,r, £) is equivalent to z being a solution of 
MiICP(G, H, F,R8,). 

(iv) Let Z = (%,u) = (x — ||ulle,w) then by using the notations and conclusion in (iii), 
we have F'(x,u) L G(x,u). We also have that 


F(x,u) = z — ||ule = z 
and 
G(x,u) = Ax + Bu + y — ||Cx + Du + vlle? 
= A(z + |lulle) + Bu + y — ||C (z + |lulle) + Du + v]le? 
= G (z, u). 
Thus, z L G(z, u). 
From the proof of (iii) we get 
0 = H(z,u) = ue! (Az + But y) + |jul|(Cz + Du + v) 
= ue! (A(z + |lulle) + Bu + y) + |lul|(C(z + lulle) + Du + v) 
= H(z, u) 
Hence, z = (x, u) being a solution of LC P(T, r, £) is equivalent to z = (x — |lulle, u) being 


a solution of MiC P(G, H, RE). 
(v) If z = (x,u) is a solution of the linear complementarity problem LC P(T, r, £) we 


have d 
z x+ Bu+y 
L> C) 1 Cae) eM. 


From (x,u) € £ and (Ax + Bu+y,Cx+Du+v) € M we have 


Tı — T2 (Ax + Buty) 
; T2 — T3 52 (Ax + Bu+y)i 
= ; € R? and : ER. 
€) Lp — Tp 522 (Ar + Bu +y); 

Lp — llull (Az + But y)i 


We also note that, from Proposition 2, it follows that for an arbitrary vector (x, u,y,v) € 
C(L), we have the following conditions 


Then, in our case, since (x, u, Ax + Bu + y, Cx + Du + v) € C(L), we have 
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%1— T2 (Ax + Buty)1 


T2 — T3 y(n + Buty); 

R4 > : =a l1 b= E€ RÅ. (6) 
at — Tp Pol (Ag + Buty); 
Lp — |lull 440+ Buy); 


where 
P 
£p = |lu], Cx + Du +v = —àu, and X (4r + Bu + y); = ||Cx + Du + v|]. 
i=1 


Then the problem of finding a solution z = (x,u) of the linear complementarity problem 
LCP(T,r, L) is converted to a problem of finding a vector z = (x,u) such that (a, 8) € 


C(R‘,). 
Moreover, let w € R? such that w; = 2; — vi41 for any i = 1,2,...,p— 1 and wp = 
Ly — ||u|| = 0. Then we have x = x(w) where x;(w) = = Wj + £p = ae w; + ||u]| for 


any i = 1,2,...,p—1 and z,(w) = |lu||. Thus, (6) is equivalent to 


Wi (Az(w) + Buty) 
we Yin (Aa(w) + Bu +y): 
R4 > : =a B= : E€ R}. (7) 
Wp-1 7-1 (Ax(w) + Bu +y): 
Wp ? (Ar(w)+ Bu +y): 


We also have from the solution of (iv) that 


H(w,u) = ue! (Ax(w) + Bu + y) + |lul|(Ca(w) + Du + v) = 0. 


Hence, the solution of (7) is equivalent to the solution of MiC P(Ĝ, Ñ, R®"'). 

(vi) Note that the function H(w, u) is a semi-smooth function and it is not differentiable 
at u = 0. Thus, we need to reformulate this function to make sure it could be differentiable 
everywhere. Let t = |lu||. Then, similarly to the proof of (v), for any (x,u, Ax + Bu + 
y, Cx + Du +v) € C(L), we have 


Tı — T2 (Ax + Buty) 
T2 — T3 yo, (Ax + Buty); 
A) =a A E€ Ri. (8) 
Tp—1 — Tp SHAS + Buty); 
Lp —t (Az Buty): 


where 


P 
Lp =t, Cx + Du +v = —\u, and X (Az + Bu +y); = |Cz + Du + v|]. 


i=1 


Next, let #2 € R? such that ù; = x; — xi+ı for any i = 1,2,...,p — 1 and ùp = £p — t = 0. 
Then, we have x = xz(ù,t), where z;(ù) = 3a Ùj + Tp = De w; +t, for any i = 


1,2,...,p — 1 and z,(w) = ||u]|. Thus, (8) is equivalent to 


w (axu, t) + Bu+y) 
w E - (Axr(ù, t) + Buty); 
R?! > k — Ig p= di (Aal q v) eRe (9) 
Wp-1 PT (Aa(ti, t) + Bu+ y)i 


We also have from the solution of (v) that 


~ ue'(Aa(w,t)+ Bu +y) +t(Cx(w,t) + Dutv 
Aout) = ( ane £2 ai 8) ) 


Hence, the solution of (9) is equivalent to the solution of MiC P(Ĝ, H, R’'). 


4 F-B function 


From the conlcusion in the Theorem 5, we have shown that the linear complementarity 
problem on the monotone extended second order cones can be converted to a Mixed com- 
plementarity problem defined on the non-negative orthant (which is defined by Facchinei 
and Pang, see Subsection 9.4.2 in [5]) is important, since by using this transformation 
scheme, the converted problem, which is the mixed complementarity problem on the non- 
negative orthant can be well studied by using the Fischer- Burmeister function, which was 
introduced by Fischer in [9,10]. For arbitrary numbers a and b, the Fischer-Burmeister 
function is defined as follows 


ola, b) = Va? + b? — (a+ b). 
From the definition of Fischer-Burmeister function, we can conclude the following property 


ola, b) = 0 <= a > 0, b > 0 and ab=0. 


We also note that ¢(a,b) is a continues differentiable function on R? \ O. By using the 
function above, for any continuously differentiable function G1, G2,...,Gp, where G = 
(Gi, G2,...,Gp), the mixed complementarity problem MiC P(G, H, R?) is equivalent to 
the following root finding problem for (x) = 0, where 


(z1, Gy) 
(x2, G2) 
B(x) = : 
(Xp, Gp) 
A 


Meanwhile, the natural merit function, 
1 2 
U(x) = 5l) 


is also continuously differentiable, and equals to zero at a point x* if and only if x* is a 
solution of MiC P(G, H,R”). Then it is equivalent to the problem of finding the stationary 
point x* of the unconstrained problem {min U(x)}. De Luca et al. used this reformulation 
to propose an algorithm that is proven to be globally convergent and locally Q-quadratically 
convergent based on considerably weaker regularity conditions than those required by the 
NE/SQP method in [4]. 

In our case, note that P, is not differentiable at u = 0, then we need to do some relaxation 


5 Generalized Newton Method for semismooth func- 
tion 


Now let Dr C R” denote the set of points at which F is differentiable. Our aim is now 
to introduce several objects from nonsmooth analysis which provide generalizations of the 
classical differentiability concept. We start by defining the B-subdifferential, where B 
stands for ” Bouligand”, who introduced the concept. 


Definition 5. Let F: U C R” > R”, where U is open and Lipschitz continues for any 
x E€ U, then the B-differential of function F at x is given by: 


OpF (x) :={GER™™ : Hrer} C Dr with xy > x, VE (x) > G} 


Algorithm 1 Newton’s method for nonsmooth systems 
1: Given F : R” > R” be locally Lipschitz continues and x, € R”, set k = 0 
2: Unless the stopping criteria is satisfied, solve the following system 


and obtain the value of dp, where G(x) is an arbitrary element of OF (xp) 
3: Set Uri = Zk + dy, k =k+1 and go back to (1). 


First, let us define the following matrix. Let Dı = diag(di1(a, u,t),..., dp—1,p—1(2, u, t)) 
and Dy = diag(d\,(x,u,t),...,d,-1-1(&, u,t)), where 


dis 2 =i 
2? + (Čr, ut) 
and ~ 
ri SH u,t) =i 


when x; Æ 0 # (Fi). Note that we also have (d; +1)? + (di, +1)? =1, then 2; = 0 = (G);, 
we have 
(da, dis) € (m2): w+ I+ (+1) = 1}. 


Then the generalised Jacobian of the FB function is the set given by 


re ol ae ee : o (10) 
Then, for an arbitrary element in the set of generalised Jacobian 
G € (x,u, t) 
when z; 40 4 G,(z, u,t), we have 
(Gz):(x, u, t) = duet + di, J_(G;(a, u, t)) 
= a -lle-4 (C — 1 | J,G(a, u,t) 


yz + (&)2(x,u, t) yz + (G)?(z, u, t) 


when x; = 0 = G;(z, u, t), we have 


(Gail, u, t) = f (dne! +d JoCla,u,t) : (du, d1) € Ball((-1,-1), 1))} 


6 Finding the minimizer of the merit function 


x is a solution to the mixed complementarity problem if it is a solution of the function 
W(x) = 0. Since U(x) is a quadratic nonnegative function, the if x is a solution of U(x) = 0 
then z is a global minimizer of function U(x). Then the problem of finding the solution to 
the mixed complementarity problem is equivalent to the problem of finding the satationary 
point of Y(x). Consider the following index sets 


C = {i: v; > 0, Hilu, v) > 0, viHilu, v) = 0}, (complementarity indices) 
R= {L2 ht C (residual indices) 

P = {i E R: v; > 0, H;(u, v) > 0}, (positive indices) 

N = {1,2,...,0Q\ (CU P), (negative indices). 


and for any arbitrary vector z, denote zs be the i-th coordinate of z, where i is an arbitrary 
number such that 7 € S and S € {C,P,N}. Then we have the following definition 


Definition 6. For the general formula of mixed complementarity problem MiC P(G, H), 
for arbitrary x € R?, u € RY and t € R, denote ŭ = (u,t)', then a point (x, u,t),is called 
FB-regular if JuG(x,u,t) is non-singular and if for any non-zero vector z € RI such that 


zc = 0, zp > 0, zy < 0, 
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there exists a non-zero vector w E€ RI such that 


and 
z' (M(z, u,t)/JeG(z, u,t)) w > 0 


where Jal ) 
xz, u, t 
G( 


) 
= J,G(2, u, t) JaG(z, u, t) ROHI) (p+a+1) 
JH (x,u,t) JaH (z, u,t) 


and M(x,u,t)/JgH(x,u,t) is the Schur complement of JzH(x,u,t) in M(z,u,t) 


In our case, for the mixed complementarity problem MiCP (Gù, u, t), A, ut), the 


Jacobian of G and H are given as 


where 
au Kis: + a12 Te Rae + ait oc igi 
sy 1 Ai j—-1 (4i + a; aa j—-1\ Gi + aig +... + Ai p- 
T dint i Peal í 2) i ini (Ga ' p-1) a eee 
i Qil DE CAI + a;2) He DET + ai +... + ent) 
where 
10. 0 
11. 0 
"a , SE 
1 1 1 
1 1 1 
01... 1 
Urs |. 2b 
0 0 1 
and A; j is a sub-matrix of A, where 
Qil Q12 Qij 
Is _ Q21 os Q2j 
Qil i2 Qij 


Pu bia avait big pu + aiz +... + aip 
m a bi = b; bgi ne bi = Qi Qi eik Qi 
B= i= 1 y= 2 > =1 “iq > =1( 1 j p) = (LıB LyAp-1p€) 
ba ET big... ORT Big SUPT (aa + Gig +... + aip) 


m * T Ax 
a= (° + ue *) 


0 
where 
Qail Qi Fai ... ss 1 aii 
A= ne a21 j Q22 «.. Xi 1 Ĉ2i 
Gp1 Qp T p2 .«-. s 1 Api 
e and 
Cii Cli T2 ~.. Ee; 1 Cli 
> c= = C21 j C22 f . i C25 
C Cgi Cgi F Cg --- 2a Cqi 
C Or equivalent to 
la ~ ~ + Ta 
> C= 
he 0 
os Moreover 
p- (P+ ue' B + (Ax(w, o + Bu +y)' elgxq Du +v + Czr(ù,t)+tCe + ue' Ae 
= = oy" 2t 
E. : where 
Le WÚ + Ú+... + wp- +t 
Ww +... + Wp- +t 
z(ů,t) = S 


Wp—ı tt 
Then, if D is non-singular, the Schur complement of D of the matrix M (W, u, y) is 


(0/5) =~ A-~BD"6 


Proposition 6. The matrix M(x,u,t) is nonsingular for any z = (x,u,t) € R? x R1 x R 
if the corresponding matriz A and D are nonsingular. 


Then we also can conclude that the Jacobian 


80(:b,u, t) = (< D) 


is non-singular if the corresponding matrices A and D are non-singular. 
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Proposition 7. Since we have the formula of Jacobian as 


jan Dı + DoJeG(w,u,t) DaJu Gld, u, t) _(D +D2A D2B 

= Ja H(w, u, t) Jun H (ô, u, t) Q D J’ 
then it is non-singular if and only if both of matriz A and D are non-singular for any 
vector z = (w,u,t) € R?*. 


Proof. For the generalised Jacobian in our case, which is 


90 (1b, u,t) = é Pe a 


We can conclude that 0®(w, u, t) is non-singular if and only if both of the Schur complement 
of D and the sub-matrix Dı + D2A are non-singular. Which is equivalent to both of the 
matrix A and D are non-singular 


The following theorem was introduced by Facchinei and Pang, for the sake of completeness, 
we quote Theorem 9.4.4 in [5] and provide a detailed proof here. 


Theorem 8. For arbitrary vectors w € R®™}, u € R? and t € R, we have z = (w,u,t) is 
a solution of MiCP(G, H, RẸ} 1) if and only if z is a FB-regular point of U(x) as well as 
a stationary point of ®(x) 


Proof. Firstly, suppose z = (w,u,t) is a solution to MiC'P(G, H, RY H. Then we have 
z = (W,u,t) is a stationary point as well as the global minimum of the associate merit 
function ®(x). Moreover, z = (Û, u,t) is a solution to MiCP(G, H,R*_'), which implies 
that (#,G(z)) € C(R%"). Then we have w = we. Thus, the FB-regularity holds for w 
and P= =N. 

Conversely, if z = (wW,u,t) is a stationary point of the merit function V(x), then 
VW(z) = 0 which implies that 


Oaou Aliut = (PH pr) Olt) =0 


Thus, for any arbitrary vector x € R?*?, we have 


= fDi +A’ Dy ČT . 
zu ae = 0. 11 
( B' D, D Pù, u, t) 0 ( ) 


For vector x we have that 


ze = 0, tp > 0, ty <0. 


Then if z is not a solution to MiCP(G, H, RY Dr we have {1,2,....p+q}\C #0. Let 
y := D(x) and we have 


yc = 0, yp > 0, yn < 0. 
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By using the definition of Dı and D2, we conclude that D,®(x) and D2®(zx) have the same 
sign. Thus, 
a'(D,®) = xé (Di®)c + xp(D1®)p + ty,(D1®)y > 0, 


since T11,9,..p+q}\c # 0, and 
a! IG(z)'(D,®)y > 0. 


Then these two inequalities above together are contradict to the condition (11). Thus, we 
have set {1,2,...,p+q}\C =9 and z is a solution to MiCP(G, Ñ, RT?) 


7 A Numerical Example 


In this section, we will give a numerical example to the liner complementarity problem 
defined on the MESOC, which is a more general case and satisfy the item (iv) in Proposition 
4. Let us consider the linear complementarity problem on the MESOC £ C R? x R?. Then 
for any arbitrary point z = (x,u) € R? x R?, the aim of finding the solution to the linear 
complementarity problem, is to find z = (x,u) € R? x R? such that (z, Tz +r) € C(L). 
By using item (vi) in Theorem 5, the solution z = (x,u) of the linear complementarity 
problem LC P(T, r, £) is equivalent to the solution of the mixed complementarity problem 
MiCP(G, H,R**) and we will have 


wy Gi (wb, u,t) (Ax(w,t) + Bu+y); 
R?! 5 Wa L Ela, u, Js (ralui, ti, t) = Xi (Arl, t) + Buty): È RP. 1 
Wp-1 Gal, u, t) PN Ar(w,t) + But y)i 
and E i . 
Ais « (Ax(w,t) + ”- v) T + Du + " 2% 
where 


iy tig t... + pitt 
thot... + tipi tt 
a(w,t) = : 
@patt 
t 


In order to finding the solution to the mixed complementarity problem, we will have the 
corresponding FB-based equation 


Let us consider the following example, where 


1 0 2 1 3 
-2 6 -1 0 -1 

t 2) = 1 -3 0 -1 -2 and r= (4) = 
0 1 -1 1 -1 
0 -1 1 1 1 


oF Fw bY 


Since we have the matrices T, A and D are non-singular, then by using the Semi-smooth 
Newton Method, the sequence {z,} = {(w,u,t),} will converge to a numerical solution to 
the mixed complementarity problem. For the solution we have 


T al 
N 82 — 12/46 „v82 — 12/46 . —225 + 30/46 139 — 24/46 
io = 0] = and u* = , . 


2 2 82 82 


Then let we check whether this solution satisfy the condition of complementarity, we have 


alg af 
82 — 12y 46 ~ 82 — 12/46 
= (EEE o) soewe (oS) 20 


Then we have 


R? > w* L Gl, u,t) € RÈ 


Then we confirm that Gi Gl, u,t)) € C(RÈŻ). 
Then by using item (vi) in Theorem 5 again, we have the the solution to the linear 
complementarity problem, which is 


i 
82 — 12V46 V82—12V46 —225 + 30V46 139 — 24/46 
“= (a= ( eee a vi v82—12/16 -225 + 30/16 m) 


2 82 82 
By using the definition of the monotone extended second order cone, we have z* € £, and 


82 — 12/46 eas 


1 0 —2 1 3 2 41 
82—12V46 2446+41 4/ 82+12V46+107 
—2 6 —l1 0 -l1 EEE 3 32 
Tz+q=| 1 -3 0 -1 -2 82—12,/46 +] 1 | = | 73v46-411/82-12V46—421 
0 1 =1 1 = 225430 /75 4 _3646ivi 
0 -l 1 1 1 139- 4V5 5 3044 6/6 
82 82 


Then by using the definition of the the dual cone of the monotone extended second order 
cone we have Tx +q E€ M and (z,Tz +q) = 0. Thus, z* is a solution to the linear 
complementarity problem. 
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8 Example for Portfolio Optimization 


As Facchinei and Pang summarized in [5], the Fischer-Burmeister function and the Gen- 
eralized Newton Method can be used to solve both the linear and the nonlinear com- 
plementarity problems. In this section, we will consider implementing this algorithm to 
solve a specific nonlinear complementarity problem, which is an application of a portfolio 
optimization problem related to the monotone extended second order cone. 

Markowitz developed the mean-variance (MV) model in [19], which is the classical 
method in investigating the problem of portfolio optimization. Suppose we build portfolio 
by using n arbitrary assets. Let w € IR” denote the weights of the assets, r € R” represent 
the return of assets and X € R” x R” be the covariance matrix. Then, the two traditional 
and equivalent MV models could be given as: 


min {wl Sw rw >a, elw= 1} 
w 


and 


max Tr'w : w! Dw < B, e'w = 1} , 
W 


where qa is the minimum profit that the investor demands and £ is the minimum risk that 
the investor wants to tolerate. They are typical quadratic optimization problems with 
higher computational complexity. 

In order to reduce the complexity of solving the portfolio optimization problem, based 
on the traditional mean-variance model, lots of models have been introduced, such as MAD 
model , which has been introduced in [15], has reduced the computational complexity 
significantly [16, 17]. 

In order to measure the uncertainty of the returns of the assets for 7 = 1,..., T, let 
us define U = (U;,...,Ur)', where U; = Ri —r. Let y; denote the upper bound of 
disturbance of return at day j. Then, the traditional MAD model can be represented as 
the following linear programming problem: 


min cf 'y—r'w 
Y,w 
T ; 
st. yj; 2 |U; wl, j=1,...,T, 
elw= 1, 


where co > 0 is the Arrow-Pratt absolute risk-aversion index. 

In reality, the uncertainty of the returns of the assets will increase with the increasing of 
the investment horizon. Thus, it is meaningful to optimize the MAD model to make it more 
in line with the real-world market behaviour. Meanwhile, by using Cauchy’s inequality, we 
also have |U} w| < ||Uj||\|w|| for any j. Then, based on the current MAD model, we obtain 
the following related problem 


min of'y—r'w 

Y,w 

s.t. YT? Urs Ss SS IU; |w|, 
e'w=1, 
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where j* = argmin, |U; w|, for j = 1,...,T. Note that the vector 


T 
YT YT-1 Yı 
3 Te , W 
Aid hil J" 


belongs to the monotone extended second order cone Lr n. Thus, the last problem is 
equivalent to the following conic optimization problem: 


min cof'y-— me 
yY, u | gt 
s.t. s (12) 


(yr, YT-1;---; Y1, u)! E LTn, 


Let us consider the KKT-conditions of the problem above. We have the Lagrange 
function as 


L(y, u) = cof'y=r' wa j -Xow — y;-1) — 01 (yı — lul) -£ ( i Ns 
Then we have 
Cofi + 62 — 01 
dL E Co f2 + ôs a 
= 4 ou m H , 
wa j Cofr-1ı + Or — Or- 
= cofr + Or 
OL r Qu 
=— + Ge 
ðu lul] 


Thus, when the condition of e'w = 1 holds, the KKT-conditions of the problem (12) can 
be converted to the following complementarity problem. 


Yı Cofi + 02 — 01 
Yo Cof2 + 03 — 02 
pal = ji l eM 
YT-1 cofr_-1+ Or — Or-1 
YT cofr + Or 
where j* = argmin, |U A w|, for j = 1,...,7 and it is a non-linear complementarity problem. 


Proposition 9. If Tal + Tat + Be £0, by using Proposition 2 and Proposition 4 ,we 
J 
have the following properties: 


(i) There exists a > 0, such that — ry + Tat + be = —Xu. 


202208.00166v1 


chinaXiv 


(ii) co¥L, fi + 207 — b1 = 


E 


T Tal 
(iti) yr = u. 


;«||, and e'w = 1, we have u # 0. Thus, item (i) and item (ii) are 
inapplicable in the problem 12 while item (iii) and item (iv) are applicable. For item (iii), 
if we have —_7— + 4" + Be = 0, which is equivalent to 


WOj*ll Tell 
er 
I Be |U; (13) 


Moreover, to make sure such 8 exists, we must have 


— iuz =... = — iun 


— OU, = 


Meanwhile, by using e! w = 


me lul (r,e) 
1=e'w=e' - i ¢ 2) 14 
AA] Iv, oe 


From (13) and (14) we have 


_ 155-16 ) 
- = (re) (15) 
Thus, 
_ lUI llU lle = r) |l 
Substitute (16) and (15) into (13) 


Thus, following the existence of 6, re any arbitrary solution (y, u) to the optimization 
problem, we must have the following conditions 


Ti — Ti+1 
Uit1 — Ui = ——_, 


(= fi ) ° 
|| U5 


and 


Moreover, from the KKT-conditions and the definition of £ we have 


Cofi + 02 — 01 
Co f2 + 03 — 02 
0= : 
cofr-1+ Or — Or-1 
cofr + Or 
Thus, 
T-1 E z 
6, = Co (2 fı fr) , when t= 1,2,...,T— 1 (18) 
—cofr, when t= T 


Then we can substitute 041 = co oS fi- fr) into (17) and we will have that the explicit 


solution of u is given as 


(Zurt yE nEn- ei e=nr) 
sny (Elar)? -n (Ear? — 0U) 
Thus, we can get the explicit solution to Problem 12. Then by using the definition of u, 


we can obtain the weight allocation of assets of the portfolio. 
Last, let we Conse = the general case. 


Suppose — rl T ye Tal “+ Be #0, then by Proposition 2, we have 
r Qu 
— be = Xu 19 
-Tel A 
and 
0u 

D -|z sie | 20 
> ‘Wit Tel oy 


i=1 
Substituting (19) into (20) we have 


T 


co X fi + 28r — 0, = Àllull. 


i=1 


Thus, 
1 T 
= Teli (a `. fi + 20r = a) 5 (21) 
i=1 


and (19) is equivalent to 
Ou 
o Be = All 


lul 
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, let the number of assets be n and we have 


Meanwhile, by using e' w = 1 and u = w||U;+ 


1 (r,e)  @,(u, e) 
l=e'w=e! a = ( : = 7 -nf | 
lUd AOI] A[U; llull 


which is equivalent to 
(r, e) A, (u, e) 


AIU; = eai 
lOs = Te Tel 


— np. (22) 


Then by (21) and (19) we get 


T 
r Qu u 
LLa EA pp 20s |\, 
[OT Tul °° Tel Qx Ta ) 


which is equivalent to 


cre (23) 
lull co Oy fi + 207 
Then we have 
n 2 
2 Boi Ti o | et r? _ 
J J i=l 
Thus, 
eae rmx (oh mj =n oe r? — U5 ||? (co ae fit 20r)’) 
— . 24 
: nll = 
By using (19) and (23) we have 
r w — Pe ( r ) ( 0; ) 
Au = — 0 2. — be = | — — ße |) {1- = 
[rl oo hae > oe? ky oer 
By using (22) and (21) we have 
(r,e) AO, (u, e) 
AllU; | = — 7 — np. 
lUl = to OS, +28, 7? 
Then we have 
be 
: ea 
u pa I[U;» || _ 0i co Di fi+20r B Be (25) 
1 (u,e) g (r,e) ` 
lU; = co Dai fi+20r—01 lU; E np 


Meanwhile, from (12), we have that (e, u) = ||U;- 


, then (25) is equivalent to 


wT Pe 
r a 0 ll jx ll — Be 
u Wl lco Xr- fi+20r 
ollUz = (re) i 
P | a -n 
sll + SE pe 1t ~ 
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Then let K = co X; fi + 207 — 61 we have 


u 7 TA — 0i — — Be 
[Ul + et at ~ 7 | 
Then 
yo Osa -PlU (14%) Wo, 
(r,e) — nB||U;« “ag” 


which is equivalent to 


_ -BUK + 41) — Or = BIU (5 2 yy 
(K + 01)((r, e) — nBU; I) (1 | | llU; 


which is equivalent to 


r= alll BN 
(K + 61) ((r, e) — nf||U; ||) (1 + | lU; 


which is equivalent to 


r — llU; le 
u= Or 
(re) — nBU 
Finally we got the weights of assets, 
we U _ F= p||U;+ e 
[U5-l| (r,e) — n BU; 


where £ is given by (24). 


9 Conclusion 


In this paper, we illustrated a method of solving a linear complementarity problem on 
the monotone extended second order cone. We have shown that the linear complemen- 
tarity problem on the monotone extended second order cone can be converted to a mixed 
complementarity problem on the non-negative orthant and reduces the complexity of the 
original problem. We can determine a solution of the mixed complementarity problem by 
using the proposition about stationary points and F-B regularity. The connection between 
the linear complementarity problem on the monotone extended second order cone, and 
the mixed complementarity problem on the non-negative orthant is also useful for appli- 
cations to portfolio optimisation. The method we illustrated works for both linear and 
non-linear complementarity problems. We expect that this scheme will also be useful for 
other applications of complementarity problems. 
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